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ABSTRACT 


The  general  linear  programming  problem  is  considered  in  which  the 
coefficients  of  the  objective  function  to  be  maximized  are  assumed  to  be 
random  variables  with  a  known  multinomial  distribution.  Three  deterministic 
reformulations  involve  maximizing  the  expected  value,  the  o'-fractile 
(or  fixed,  0  <  a  <  ^),  and  the  probability  of  exceeding  a  predetermined 
level  of  payoff,  respectively.  In  this  paper  the  author's  previous 
work  on  "bi -criterion  programs"  is  applied  to  derive  an  algorithm  for 
routinely  and  efficiently  solving  the  second  and  third  reformulations.  A 
by-product  of  the  calculations  in  each  case  is  the  tradeoff-curve  between 
the  criterion  being  maximized  and  expected  payoff.  The  intimate  relationships 
between  all  three  reformulations  are  illuminated. 


ON  STOCHASTIC  LINEAR  PROGRAMMING 


A„M.  Geoffrion 

1.  INTRODUCTION 

Consider  the  problem 

(1)  Maximize  p^x  subject  to  Ax  £  b  , 

vhere  x  is  an  n-vector,  p  is  a  random  n- vector  with  a  known  distribution, 
and  the  linear  constraints  are  deterministic.  x  must  be  chosen  knowing 
only  the  distribution  of  p.  Three  deterministic  reformulations  of  (l)  are: 

(2)  Maximize  E(x)  subject  to  Ax  <  b, 
where  E  denotes  the  expected  value  of  p^x; 

(3)  Maximize  F  (x)  subject  to  Ax  <  b, 

a  — 

where  0  <  o  <  1  is  a  predetermined  constant  and  F  (x)  is  the 

a 

~t 

of-fractile  of  p  x;  and 

(*0  Maximize  P^(x)  subject  to  Ax  <  b  , 

where  P^(x)  is  probability  that  p^x  equals  or  exceeds  a  predetermined 
"aspiration"  level  h  of  payoff.  The  expected  value  reformulation  (2)  has 
the  computational  advantage  that  it  leads  to  an  ordinary  linear  program  to 
be  solved,  whereas  the  fractile  and  "aspiration"  reformulations  —  which  can 
be  more  realistic  in  certain  situations  —  lead  to  nonlinear  programs. 

The  purpose  of  this  paper  is  two-fold:  to  apply  the  author’s  previous 
work  cn  "bi-criterion  programs"  go  as  to  obtain  an  algorithm  for 
routinely  solving  (3)  or  (4)  by  parametric  quadratic  programming;  and  to 
illuminate  the  intimate  relationships  between  (2),  (3)  and  (4),  The  algorithm 
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presented  here;  although  turning  out  to  have  much  in  common  with  certain 
previous  approaches  to  (3)  or  (4),  offers  a  unified,  rigorous,  non-graphical, 
and  computationally  efficient  approach  to  (3)  and  (4).  It  has  the  additional 
advantage  of  yielding  as  a  by-product  of  the  calculations  the  tradeoff  curve 
between  the  criterion  function  being  maximized  and  E(x). 

For  simplicity  we  assume  that  p  is  multinormal  with  mean  vector 
p,  and  covariance  matrix  2,  so  that  p^x  is  N(pSc,  x^Ex);  that  (2)  has  a 
finite  optimal  value  Mj  and  that  x^Ex  >  0  for  all  feasible  x. 

We  conclude  this  Introduction  with  some  preliminary  remarks  on  the 
tractile  and  aspiration  criteria,  including  a  survey  of  known  computational 
approaches. 

The  Practile  Criterion 

Since  the  .5-fractile  (i.e.,  the  median) of  a  normally  distributed 
random  variable  coincides  with  its  mean,  (3)  with  a  =  .5  is  identical  with 
(2).  Maximizing  the  a-fractile  with  0/  <  .5  should  appeal  to  some  conserva¬ 
tive  decision-makers  because  it  tends  to  control  the  lower  tail  of  the 
distribution  of  payoffs.^/ 

It  i3  easy  to  show  that 

Fa(x)  =  nSc  +  /x^Ex  , 

where  $  is  the  Standardized  Normal  Distribution  Function.  When 
0  <  a  <  .5,  as  we  assume  henceforth,  <  0  and  (3)  is  a  concave 

program,  since  /x^Ex  is  convex  [9;  P.  195].  Note,  however,  that  the 
criterion  function  is  not  differentiable  when  x^Ex  vanishes,  thus  limiting 

1 /  Cf.  Baumol  [33,  who  seems  to  be  getting  at  this  idea  in  the  context 
of  the  portfolio  selection  problem,  which  is  a  special  but  important 
case  of  (l). 
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the  applicability  of  gradient-type  optimization  procedures.  Eisenberg  [73* 
Sinha  [l2j,  and  others  have  stepped  into  the  breach  with  theoretical  results 
that  are  designed  to  facilitate  a  computational  solution. 

Computational  procedures  for  (3)  have  been  offered  by  Kataoka  [9], 
who  proposed  and  partially  justified  an  iterative  procedure  that  can  be  viewed 
as  a  discretized  variant  of  the  algorithm  derived  here,  and  by  Sinha  tl2), 
who  developed  an  elaborate  specialized  duality  theory  that  leads  to  a 
computational  solution  involving  linear  and  quadratic  programming.  It  can 
also  be  solved  by  an  obvious  geometrical  construction  that  requires  the  graph 
of  the  (E,<j)-  tradeoff  curve  arsociated  with  (l),  i.e.,  the  image  of  all 
feasible  x  with  the  property  that  a  higher  value  of  can  be  attained 

only  at  the  expense  of  a  higher  value  of  /“x^ST,  and  a  lower  value  of 
/"xt-ISc  only  at  the  expense  of  a  lower  value  of  jiSc.  The  (E,o) -tradeoff 
curve  is  most  conveniently  obtained  by  a  square-root  transformation  from 
the  (E,ct  ) -tradeoff  curve.  The  latter  should  be  computed,  as  pointed  out 
by  Markowitz  [lO],  by  parametric  quadratic  programming.  This  use  of  parametric 
quadratic  programming  is  at  the  heart  of  the  present  approach  to  (3)  and 
(4),  but  we  do  not  require  a  graphical  construction  and  only  compute  a 
relevant  subset  of  the  (E,ct^)-  tradeoff  curve. 

The  Aspiration  Criterion 
Clearly 

pfx)  -  1  -  *(  22!*-)  s  (  (  j±^=a_)  . 

/  xt£x  /  x^Ex 

Since  $  is  strictly  increasing,  (4)  has  the  same  optimal  solution  set  as 


(*A) 


M^x-K 


subject  to  Ax  <  b. 


Maximize 


/  xt2x 
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When  the  aspiration  level  k  is  taken  as  M  (the  optimal  value  of  (2)), 
it  is  easily  seen  that  (4a)  has  an  optimal  value  0  and  that  this  value  is 
achieved  for  those  feasible  x  for  which  m.^x  =  M.  Hence  (4a)  with  h  =  M 
has  the  same  optimal  solution  set  as  (2),  and  the  maximum  probability  in  (4) 
is  (since  $(0)  =  i).  As  before,  we  are  interested  in  the  conservative 
decision-maker,  who  would  undoubtedly  take  k  <  M.  Such  a  choice  is  necessary 
and  sufficient  for  the  maximum  probability  in  (4)  to  be  >  0.5>  and  we  assume 
it  henceforth. 

The  aspiration  criterion  program  (4a)  can  be  solved  by  a  simple  geometrical 
construction  noticed  by  Roy  Cll],  who  presented  it  for  a  special  case,  providing 
that  the  (E,cr)-  tradeoff  curve  is  available.  The  method  can  be  modified  to 
work  almost  as  easily  with  the  (E,a2)-  tradeoff  curve,  which  as  we  pointed 
out  above  is  easier  to  compute. 

See  Charnes  and  Cooper  [4]  for  additional  discussion  of  the  aspiration 
criterion.  They  also  show  how  to  reduce  (4a)  to  a  program  that  is  linear 
except  for  one  quadratic  constraint.  If  this  quadratic  constraint  were  dealt 
with  by  the  standard  trick  of  taking  it  up  into  the  objective  function  with  an 
undetermined  multiplier  and  applying  quadratic  programming,  the  result  could 
be  an  algorithm  that  closely  resembles  the  one  given  here. 

2.  THE  ALGORITHM 

In  this  section  we  view  (3)  and  (4)  as  if  they  were  bi-criterion  programs 
and  derive  an  algorithm  for  each  by  applying  the  following  result  from  [8], 

Let  it  be  desired  to  solve 

(5)  Maximize  ufp^f^x)),  P2(f2(x))  ), 
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where  X  is  a  non-empty  compact  convex  set  in  Rn,  f-^,  f2>  p2(f2^ 

are  concave  on  X,  p^  and  p2  are  strictly  increasing  on  the  image  of 
X  under  f^  and  f 2,  respectively,  u  is  non-decreasing  and  quasiconcave 
(see,  e.g.,  [l])  on  the  convex  hull  of  the  image  of  X  under  (p^( f^),  P2(f2))# 
and  all  functions  are  continuous.  Think  of  u(»,*)  as  a  utility  function 
defining  a  preference  ordering  over  pairs  of  values  of  the  two  criterion 
functions  f^  and  f 2,  on  which  the  scale  transformations  p1  and  p2  have 
been  performed.  Assume  that  a  parametric  programming  algorithm  is  available 
for  solving 

(6)  Maximize  y  f.(x)  +  (l-y)  f?(x) 

xeX  1 

for  each  value  of  the  parameter  y  in  the  unit  interval,  and  that  the 
resulting  optimal  solution  function  x*(y)  would  be  continuous  on  [0,1]. 

Then  the  function 

u(y)  *  u(px(  fx(x*(  v))  ),  p2(f2(x*(v))  )  ) 

is  continuous  and  unimodal  on  Co,  1 3,  and  if  y*  maximizes  U(y)  on  [0,1] 
then  x*(y*)  is  optimal  in  (5). 

Solving  (3) 

Consider  now  (3).  Put 
X  =  (x:  Ax  <  b} 
fx(x)  =  p^x 
f2(x)  =  -xtSx 

Pl(fl>  =  fl 
•^2^  ^2^  =  ~  ^ 

u(pl(fl)>  P2(f2))  =  P^^)  “  P2(f2^ 
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It  is  easy  to  verify  that  all  of  the  assumptions  required  of  (5),  which  in 

this  case  is  identical  to  (3),  are  satisfied  save  one:  the  compactness  of  X. 

certain 

However,  the  only  need  for  compactness  is  to  ensure  that  /  supreme  are 
achieved.  The  attainment  of  all  supreme  here  follows  from  the  non-negativity 
of  /~&Tx.  ,  our  assumption  that  (2)  has  a  finite  optimal  value  M,  and 
the  fact  [2,  Th.  1.73  that  a  concave  quadratic  polynomial  bounded  above  on 
a  convex  polyhedral  set  achieves  its  constrained  supremum.  Now  (6)  becomes 
(7)  Maximize  y  (ji^x  -  (l-y)  x^Ex  subject  to  Ax  <  b, 

a  parametric  quadratic  program.  Several  algorithms  are  available  for  it 
(e.g.,  [lO],  [l3],  [5])#  and  they  all  yield  an  optimal  solution  function 
x*(y)  that  is  continuous  on  [o,l]«  Hence  the  method  quoted  above  applies, 
and  we  see  how  to  solve  (3)  with  the  aid  of  any  parametric  quadratic  programming 
code  for  (7). 

To  solve  (3)  one  may  solve  (7)  with  y  =  and  decrease  y  until  the 

unimodal  function  p  (x*(y))  achieves  its  maximum  on  [0,l].  When  the  maxi- 

Qf 

mizing  y  is  reached,-^  say  Ya  ,  the  parametric  programming  is  stopped 
because  the  optimal  solution  x*(  y^ )  of  (3)  has  been  found.  For  more  complete 
details  see  [8],  where  the  fact  that  the  image  of  lo,l]  under  x*(y)  is 
piecewise  linear  is  taken  advantage  of  in  an  obvious  way. 

Solving  (4) 

The  assumption  h  <  M  guarantees  that  the  optimal  value  of  (^A)  is  >  0. 
Therefore  one  can  restrict  attention  in  (4a)  to  feasible  x  such  that 

2/  For  a  reason  that  will  become  apparent  in  the  next  section,  y  =  1  is  a 
natural  storting  point  --  although  any  other  value  in  the  unit  interval 
could  be  used. 

3/  If  the  maximizing  value  is  not  unique,  let  ^  be  the  largest. 
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IJ^x  -K  >  0.  In  thi5  region  the  maximand  is  quasi -concave.  Put 

X  =  {x:  Ax  <  b,  p/tx-H  >  0) 

f^x)  =  p^x 
f2(x)  =  ^Ix 

pl(fl)  =  fl"H 

p2(  =  "  v/Cp2 

uCp^^),  p2(f2))  =  px( fx)/  (-P2(f2)). 

Ic  is  easy  to  verify  that  all  assumptions  required  of  (5)*  which  in  this 
case  is  identical  to  (4a),  are  satisfied  save  the  compactness  X.  As  before, 
this  seeming  difficulty  is  eliminated  by  the  fact  that  all  pertinent  constrained 
suprema  are  achieved.  Furthermore,  (6;  is  again  the  parametric  quadratic 
program  (7).“-'  Thus  to  solve  (4a)  one  may  solve  (T)  for  Y  =  1  and  decrease 
V  until  the  unimodal  function  (p."tx*(  y)-x)/  /x*(  y)£x*(  y)  reaches  its 
maximum  on  [o,  1  ] .  When  the  maximizing  y  is  reached,-^  say  Y  t  the  optimal 

H 

solution  x*(y  )  has  been  found.  Again,  consult  [8]  for  details, 

H 

3.  INTERPRETING  THE  INTERMEDIATE  QUANTITIES 
The  algorithm  for  solving  (3)  or  (4)  given  in  the  previous  section 
involves  the  computation  of  x*(y)  (an  optimal  solution  of  the  parametric 
quadratic  program  (7))  from  right  to  left  on  the  interval  # l]  in  one  case 
and  on  Lyh>  l]  in  the  other.  Although  x*(Ya)  ^  x*(Yk)  63:6  °Pti®al  solutions 

hj  Technically,  (7)  must  now  include  the  constraint  JjI^x-h  >  0.  Since  we 
shall  take  y  to  be  decreasing  from  1,  this  constraimTwill  never  be 
binding  before  termination  and  can  therefore  be  dropped. 

5/  If  the  maximizing  value  is  not  unique,  let  y„  be  the  largest. 
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to  (3)  and  (4),  the  intermediate  x*(y)  are  also  of  interest.  We  shall  show 
in  Theorems  3  and  4  of  this  section  that  the  image  in  R  of  [y^jI]  under 
(e(x*(y))>  Fa(x*(Y))  )  is  the  complete  (E,  F^)-  tradeoff  curve,  and  that 
the  image  in  R2  of  [y  ,lj  under  (E(x#(y))*F.  (x*(y))  )  is  the  complete 
(E,  P  )-  tradeoff  curve. We  also  show  in  Th.  2  that  each  intermediate 

H 

x*(y)  also  solves  (3)  with  a  certain  a  and  (4)  with  a  certain  h. 

It  will  be  notationally  convenient  in  the  sequel  to  refer  to  the  following 
reparameterized  version  of  (T): 

(7A)  Maximize  p,^x  -  Sx^Ex  subject  to  Ax  <  b, 

where  the  parameter  3  traverses  [0,°°).  For  fixed  3  e  to,00),  (7A)  is 
equivalent  to  (7)  with  Y  =  (1+8)”1,  and  therefore  has  an  optimal  solution 
x*(  )  which  with  some  abuse  of  notation  we  henceforth  call  simply  x#(8). 

Let  3a  -  (i-Va)/Ya  and  8k  =  (1-Yk)/yk.  Then  (3)  or  (4)  are  solved  by 

computing  x*(3)  from  left  to  right  on  Co,  3  ]  or  [0,8  ],  respectively. 

Of  H 

Theorem  1: 

(i) .  Let  ot  be  fixed  in  the  unit  interval.  Every  optimal  solution  of 

(3)  is  also  an  optimal  solution  of  (4)  with  h  equal  to  the 
optimal  value  of  (3). 

(ii) .  Let  k  be  fixed  arbitrarily.  Every  optimal  solution  of  (4)  is 

also  an  optimal  solution  of  (3)  with  (l-a)  equal  to  the  optimal 
value  of  (4). 

2 

6/  Naturally  x*(y)  also  yields  a  portion  of  the  (E,ct  )-tradeoff  curve,  but 
this  fact  is  of  little  more  than  passing  interest  here. 
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Proof:  Let  x  and  x*  be  feasible.  Then 

(i^x  +  /x^Hx  <  pSc*  +  §~\&)  / x*^v* 

if  and  only  if 

-  u-Sc*  -  $“1(a)  /x^Sx*  <  -  $-1(a)  /jT-Sc 
if  and  only  if 


li^x-n^x*-?  ^(a)  /x*1'!** 

/?5 


■trr 


if  and  only  if 


p.tx-((j.tx*t?  ^(ar)  <  m-^x*  -  (^x*t$”^(cv)  /x*^Xoc*) 

/x^Ex* 

Hence  for  fixed  x*  the  first  inequality  holds  for  all  feasible  x  if  and 
only  if  the  last  holds.  The  theorem  follows. 

This  theorem  points  up  a  reciprocity  between  the  fractile  and  aspiration 
criteria  that  holds  in  much  more  general  circumstances  than  the  one  con¬ 
sidered  here. 

Theorem  2: 

For  3  >  0,  x#(3)  (an  optimal  solution  of  (TA))  is  also 

(i) .  optimal  in  (3)  with  ot  =  $  (-23  /x*( 3)  ) 

(ii) .  optimal  in  (4)  with  k  =  p^x*^)  -  23x*(<3)tEx*(3). 

Proof:  Part  (i),  which  is  essentially  equivalent  to  Th.  3  in  [9],  can  be 
proven  by  setting  up  the  appropriate  identification  between  the  Kuhn-Tucker 
conditions  for  (TA)  and  for  (3).  Part  (ii)  follows  from  part  (i)  and 


Theorem  l.i 


11 


These  two  theorems  give  an  interesting  interpretation  to  the  intermediate 

x*(0).  In  the  course  of  solving  (3)  with  a  =  afQ(0  <  otQ  <  by  the  method 

of  section  2,  one  automatically  solves  ( 3)  with  each  value  of  of  between 

<*o  and  ^  and  (4)  with  each  value  of  k  between  F  (x*(0a  ))  (which 

o  o 

must  be  <  M)  and  M.  Similarly,  in  the  course  of  solving  (4)  with  k  =  k0 

(h0  <  M)  by  the  method  of  section  2,  one  automatically  solves  (4)  with 

each  value  of  k  between  h0  and  M  and  (3)  with  each  value  of  Of  between 

1-P  (x*(P  ))  (which  must  be  <  b)  and 

Ko  Ho 

Lemma: 

The  function  2fVx*(0)  ^Lx*( 3 )  is  non -decreasing  onto,  3^  ]•  It  assv 
the  value  0  only  at  0=0  and  the  value  1(ar)  only  at  (3=0^. 

Proof:  The  proof  requires  two  applications  of  a  method  often  used  for 
obtaining  similar  mcnotonicity  results  in  parametric  programming. 

Let  0^  and  B  satisfy  0  <  6^  <  3  <  . 

First  we  derive  the  preliminary  inequality 

x*(0  )t£x*(0  )  <  x*(0°)tLx*(0°). 

By  the  definitions  of  x*(0°)  and  x*(0*)>  we  have 

^*(3’)  -  0°x*(0' )*&<*(  B')  <  M>tx*( 0° )  -  0°x*(0°)t2x*(0°) 

and 

jj,^x*( B° )  -  3'x*( 0°)tlx*( 0°)  <  ^*(0')  -  0’x*(0»  )*£**(  p'). 

Summing  and  rearranging  these  inequalities,  we  find 

(0'-0°)(x*(0')t  2x*(0')  -  x*(0°)tlbc*(00))  <  0, 
which  upon  division  by  (0'-0°)  yields  the  preliminary  inequality. 
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By  Th.  2.i,  x#(P°)  maximizes  p^x  -  2p°/x#(  Bu ) ^ Ex*(  Pu )  /x^Ex  over 

all  feasible,  x.  A  similar  assertion  holds  for  x*(8').  Applying  the 
above  argument  to  such  programs  in  lieu  of  (7A),  one  obtains 

(b*  -8°v/x*r^°)tS*(80)  )  (/7»(¥7Tt^*^1TT 

-  /x*TF°)tac*(?0)  <  o. 

The  preliminary  inequality  implies  that  the  second  factor  is  either  negative 

or  zero.  In  the  first  case#  division  by  it  yields  that  the  first  factor  is 

•  o 

>  0;  in  the  second  case,  the  first  factor  is  >  0  because  3  -0  >  0.  Thus 
the  first  factor  must  be  non-negative,  and  the  desired  monotonicity  of 
2fVx*(0)tEx*(B)  is  established. 

That  the  value  0  is  assumed  only  at  6=0  follows  from  our  assumption 
that  x^Ex  >  0  for  all  feasible  x  ;  that  the  value  is  assumed 

only  at  0  follows  from  the  nature  of  our  algorithm  for  solving  (3). 

a 

The  proof  is  complete. 

Theorem  3: 

traverses  .  o 

As  0/  [0,  0^],  (p,  x*(0),  Fr/(x*(3)))  traces  in  R  the  complete 

(E,Fa) -tradeoff  curve. 

Froof:  Since  x*(0)  is  optimal  in  (2),  x*(B^)  is  optimal  in  (3),  and 

(p.Sc^KB),  Fa(x#(0)))  on  [D,  P^]  determines  a  curve  by  the  continuity  of 
x*(0),  it  is  sufficient  to  show  that  x#(0)  is  (E,  F^) -efficient  on 

t 

(0,0^).  To  show  this  it  iA  clfeprly  sufficient  to  show  that  for  0  <  0  <  Ba 
there  exists  a  scalar  X,0  <  X  <  1,  such  that  x#(3)  maximizes 
(l-X)ix^xtX  Ffl(x)  over  all  feasible  x  .  This  maximand  simplifies  to 
p.txtX$”1(a)  /xfcEx”’,  For  X=20/x*( B )^Ex*( 3 )/( -$  ^(a)),  x*(B)  indeed  max¬ 

imizes  it  by  Th.  2.1  .  The  Lemma  implies  that  this  choice  of  X  satisfies 
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0  <  X  <  1  .  Hence  x*(0)  is  (E,F^) -efficient  for  0  <  B  <  S^,  and 
the  proof  is  complete. 

Theorem  4: 


As  0  traverses  [o,  8  3,  (m^x^B),  P  (x*(P)))  traces  in  R  the 


H 


complete  (E,  P  ) -tradeoff  curve. 


H 


Proof:  Since  x*(o)  is  optimal  in  (2),  x*(P^)  is  optimal  in  (4),  and 

(^x^B),  P  (x*(0)))  on  [o,  B  ]  determines  a  curve  by  the  continuity  of 

*  K 

x*(8),  it  is  enough  to  show  that  x*(G)  is  (E,  P  ^-efficient  on  (0,  B  ).  To 

H  H 

show  thi3  we  shall  resort  to  the  definition  of  (L,  P  )-efficiency  and  demonstrate: 
for  0  <  0<  0  »  one  has 

H 

(8)  P  (x)  <  P  (x*(0)) 
for  all  feasible  x  such  that 

(9)  P-^x  >  |i^x*(8)  , 

with  equality  in  (8)  impling  equality  in  (9)  , 

Clearly  one  may  use  (^x-x)/  /x^Zx  in  place  of  P  (x)  in  (8).  Hence 

H 

(8)  can  be  written 

(8A)  (ii^x-k)/  /x^Zx  <  (hV([3)-h)/  /x*('3)tZx*(0)  . 

The  method  of  demonstration  will  use  the  theory  of  linear  programming,  after 
some  manipulation. 

Let  8  satisfy  0  <  0  <  0  .  By  Th.  2.i, 

(10)  ix^x-20  /x*(3)tZx*(P)  /x^Zx  <  \2^x*(  B)-2px*(  B)^2x*(B) 


for  all  feasible  x.  It  is  therefore  sufficient  to  show  (8a)  for  all 
pairs  (m-^x,  /x^Zx)  satisfying  (9)  and  (10),  with  equality  in  (8A) 
implying  equality  in  (9).  It  is  now  convenient  to  change  notation. 

Let  y^  =  p^x  -  u^x^B),  y2  =  /x^Zx  .  Then  the  proposition  to  be  shown  can 
be  written  as  follows  after  rearranging  (8A): 
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(8B)  y1  /x*{  3)^*0)  -  ^-(h^OJ-h)  /x*(0)  lEx*{  P)  . 

for  all  (y  ,  yg)  >  0  satisfying 

( 10 A)  yx  -  20  /x*( 0)^(3)  y2  <  -20x*( 3)tSx*(P)  , 

with  equality  in  (8B)  implying  y^=0.  The  identity  of  a  simple  linear 
program  now  becomes  apparent,  with  the  left-hand  side  of  (8b)  serving  a s 
a  linear  objective  function  to  be  maximized  and  (lOA)  serving  as  a  linear 
constraint.  This  linear  program  is  feasible  (e.g.,  take  y^O  and 
y^  =  /x*(  *-(b)  ).  The  object  then  becomes  to  show  that  the  optimal 

value  of  this  linear  program  is  equal  to  the  right-hand  side  of  (8b),  and 
that  y-^  necessarily  vanishes  in  any  feasible  (y^,y2)  pair  that  achieves 
this  value.  By  the  Dual  Theorem  and  Complementary  Sir"’  -ss  Theorem 
of  linear  programming  [6],  it  is  sufficient  to  demonstrate  tnat  the  dual 
linear  program  is  feasible,  that  its  optimal  value  equals  the  right-hand 
side  of  (8B),  and  that  the  dual  constraint  corresponding  to  y1  necessarily 
is  satisfied  with  strict  inequality  at  every  optimal  dual  solution.  The 
dual  linear  program  is: 

Minimize  (-20xJ*(0)'*'2x*(0))  z 

subject  to  z  >  /  x*(p)^£x*(c) 

(-20  /x*(P)t'Zx*CP7  )  z  >  -  (^*(0)  -  h) 


z  >  0 
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From  our  assumption  that  n  <  M,  the  continuity  of  x*(3)  on  [0,  3^],  the 
nature  of  our  algorithm  for  (4),  and  Th.  2.ii,  it  follows  that 

h  <  -  23x*(  P^Ex*^  3)  on  [o,  3  ).  With  this  inequality  it  is 

H 

now  easy  to  verify  that  z*  =  (n^x^P)  -  h)/  23  / y.*{  8)’tILx*(  S)  is  the 
unique  optimal  solution,  and  consequently  that  the  dual  program  has  the 
requisite  properties.  Thi3  completes  the  demonstration. 

Theorems  3  end  4  reveal  why  the  initial  value  of  3  was  chosen  to 
be  0.  Any  non-negative  initial  value  could  be  chosen  and  the  algorithm 
would  still  solve  (3)  (resp.  (4)),  perhaps  even  with  less  calculation;  but 
if  the  initial  value  exceeds  3  (resp.  3  ),  then  the  intermediate 

0!  K 

x*(3)  will  no  longer  be  (E,F  )  (resp.  (E,P  ))  -  efficient,. 

Ur  K 
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